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FIG. 3 
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Determine a first gas force waveform 
relating to a gas force acting on 
a first cylinder piston 




Determine a second gas force waveform 
relating to a gas force acting on 
a second cylinder piston 
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Determine a first inertia force 
waveform relating to an inertia 
force acting on a first cylinder 

piston 




Determine a second inertia force 
waveform relating to an inertia 
force acting on a second cylinder 

piston 
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Determine a first forcing waveform 
indicative of crankshaft torque 
at a first crankshaft throw, 
based on the first gas force 
waveform and the first 
inertia force waveform 




Determine a second forcing waveform 
indicative of crankshaft torque 
at a second crankshaft throw, 
based on the second gas force 
waveform and the second 
inertia force waveform 
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FIG. 4 




! 

Remove a first DC component from the first 
forcing waveform to obtain a first 
modified forcing waveform 
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Store the first DC component 
in a memory 
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Apply a Fourier transform to 
the first modified forcing waveform 
to obtain a first plurality of complex spectral 

values 




Frequency 1 1.2 +j 1.6 
Frequency 2 1.0 + j 1.5 
Frequency 3 1.5 + j2.1 
Frequency 4 0.8 + j0.2 
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FIG. 5 
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Calculate a first plurality of spectral amplitude values based 
on the first plurality of complex spectral values 

Frequency 1 2.0 

Frequency 2 1 .8 

Frequency 3 2.6 

Frequency 4 0.8 
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Determine a first maximum spectral amplitude from the first 

plurality of spectral amplitude values 

[Frequency 3 2.6] 




Determine an overall maximum spectral amplitude by calculating 
the greater of the first maximum spectral amplitude and the 

second maximum spectral amplitude 
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Determine a threshold amplitude value based on the overall 
maximum spectral amplitude and an acceptance value: 
threshold amplitude value = overall maximum spectral amplitude 

* acceptance value 
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FIG. 6 
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Copy portions of the second 
forcing waveform to itself to obtain 
an integral number of revolutions 
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Remove a second DC component from the 
second forcing waveform to obtain a second 
modified forcing waveform 
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Store the second DC component 
in a memory 
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Apply a Fourier transform to 
the second modified forcing waveform 
to obtain a second plurality of complex 

spectral values 

Frequency 5 1.2+j1.6 

Frequency 6 1.0 + j 1.5 

Frequency 7 1.5+j2.1 

Frequency 8 0.8 + j0.2 
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FIG. 7 
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Calculate a second plurality of spectral amplitude values based on 
the second plurality of complex spectral values 



Frequency 5 
Frequency 6 
Frequency 7 
Frequency 8 



2.0 
1.8 
2.5 
0.8 
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Determine a second maximum spectral amplitude from the second 

plurality of spectral amplitude values 

[Frequency 3 2.5] 
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FIG. 8 
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Determine a first plurality of desired frequency values by selecting 
frequency values associated with a subset of the first plurality 
of spectral amplitude values that are greater than or equal to 

the threshold amplitude value 



Frequency 1 
Frequency 2 
Frequency 3 
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Determine a second plurality of desired frequency values by selecting 
frequency values associated with a subset of the second plurality 
of spectral amplitude values that are greater than or equal to 

the threshold amplitude value 

Frequency 5~* 



Frequency 6 
Frequency 7 
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Solve a dynamic model using: (i) a subset of the first plurality of 
complex spectral values associated with the first plurality of 
desired frequency values, and (ii) a subset of the second plurality of 
complex spectral values associated with the second plurality of desired 
frequency values, and (iii) the first and second DC components 
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Frequency 1 1.2+j1.6 
Frequency 2 1.0+j1.5 
Frequency 3 1.5+j2.1 



Frequency 4 1.2+j1.6 
Frequency 5 1.0+j1.5 
Frequency 6 1.5+j2.1 



First and second DC components! 



End y 
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